Astronomy & Astrophysics manuscript no. 14476'cl 


©ESO 2010 


June 11,2010 





Kinematics and line strength indices in the halos of the Coma 
Brightest Cluster Galaxies NGC 4874 and NGC 4889 * 

L. Coccato 1 , M. Arnaboldi 2 ' 3 , O. Gerhard 1 , K. C. Freeman , G. Ventimiglia 1 ' 2 , and N. Yasuda 5 

1 Max-Plank-Institut fur Extraterrestrische Physik, Giessenbachstra/?e 1, D-85741 Garching bei Miinchen, Germany. 

2 European Southern Observatory, Karl-Schwarzschild-Stra/?e 2, 85748 Garching bei Miinchen, Germany. 

3 INAF, Osservatorio Astronomico di Pino Torinese, 1-10025 Pino Torinese, Italy. 

4 Research School of Astronomy & Astrophysics, ANU, Canberra, Australia. 

5 Institute for the Physics and Mathematics of the Universe, University of Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8568, 
Japan. 

June 11, 2010 

ABSTRACT 

Aims. We investigate the stellar kinematics and line strength indices in the outer halos of brightest cluster galaxies (BCGs) in the 
Coma cluster to obtain the outer halo V ro t and <r profiles and to derive constraints on the formation history of these objects. 
Methods. We analyzed absorption lines in deep, medium-resolution, long-slit spectra in the wavelength range ~ 4500 - 5900 A, out 
to ~ 50 kpc for NGC 4874 and ~ 65 kpc for NGC 4889, probing regions with a surface brightness down to ji R ~ 24 mag arcsec~ 2 . 
Results. These data provide stellar velocity and velocity dispersion profiles along the major axes of both BCGs, and also along the 
minor axis of NGC 4889. The kinematic properties of NGC 4874 and NGC 4889 halos extend the previous relations of early-type 
galaxy halos to bright luminosities and indicate that the stars in the outer regions are still bound to these galaxies. For NGC 4889 we 
also determine H/3 Mg and Fe line strength indices, finding strong radial gradients for Mg and Fe. The current dataset for NGC 4889 
is one of the most extended in radius, including both stellar kinematics and line strength index measurements. 

Key words, galaxies: cluster: general - galaxies: cluster: individual: Coma cluster - galaxies: individual: NGC 4874 and NGC 4889 
- galaxies: kinematics and dynamics - Galaxies: abundances 



1. Introduction 

Brightest cluster galaxies (BCGs) are large and luminous gala- 
xies located in the centers of galaxy clusters. The formation his- 
tory of BCGs an d their halos is connected to the formation of 
the cluster itself (Dubinski 1998) a nd to t he presence of diffus e 
intra-cluster light (iGonzalez et all 120051: iMurante etail 120071) . 
The formation of BCG and their halos can be investigated by 
the combined study of the stellar kinematics and the population 
content of their halos. Dynamical timescales in the halos are on 
the order of 1 Gyr Q, approaching a significant fraction of the 
age of the universe, and therefore the fingerprints of the forma- 
tion processes may still be preserved there. 

So far, the studies that provide both stellar kinematics and 
line strength indices in BCGs have been limited to within 
one effective radius (e.g | Carter et ail 1999t Fisher et al. 
Brou gh et al.ll2007l: ISpolaor et al. 2008a b; Loubser et al. 



1995 



2008 



2009), therefore measurements over a wider radial range are 
highly desirable. 

This work is the first of a series aimed at studying the for- 
mation history of BCGs by probing the stellar kinematics and 
populations of their outer halos, covering regions at 3 effective 
radii or larger. As first targets, we selected NGC 4874 and NGC 
4889, the two BCGs in the Coma cluster (Abell 1656). The in- 



Send offprint requests to: L. Coccato, e-mail: lcoccato@mpe.mpg.de 
* Based on data collected at Subaru Telescope, which is operated by 

the National Astronomical Observatory of Japan. 
1 Estimat ed for distance R = 50 k pc and circular velocity V c = 250 

km s _I , e.g. Binnev & Tremaine 1 9871 . 



ner parts of these galaxies have been studied with photometric, 
kinematic and stellar populations data documented and av ailable 
in the hterature (e.g. iMehlert et"al] l2000t iGavazzi etail 12003: 
Iderhard et ai1l2007r.lTrager et al.ll2008h . 

In this paper we describe the data acquisition, reduction, and 
the measurements of the long slit stellar kinematics for both ga- 
laxies, and the line strength indices for NGC 4889. The data set 
in NGC 4889 extends out to 65 kpc ( which correspond to ~ 4.3 
effective radii, Jorgensen et alj[T995l) . providing the most radi- 
ally extended measurements of absorption line kinematics and 
line strength indices in the outer halo of a BCG. These data are 
the basis of forthcoming papers investigating the formation his- 
tory and evolution of the galaxies in the Coma cluster core. 

The paper is organized as follows. Spectroscopic observa- 
tions and data reduction are discussed in Section[2] The sky sub- 
traction and the radial binning of the spectra are described in 
Sections[3]and|4] respectively. Section [5] describes the measure- 
ments of the stellar kinematics, while Section [6] describes the 
measurements of the line strength indices. Finally, the results 
are discussed in Section[7] 

In this paper, we adopt a distance to NGC 4784 of D = 102.6 
Mpc and an effective radius of R e = 70'.'79 = 35.21 kpc; for 
NGC 4889 we adopt D = 92.7 Mpc and R e = 33788 = 15.23 
kpc. Distances are taken from the NASA/IPA C Extragalactic 
Datab ase (NED) and the effective radii from Jorgense n et alj 
(Il995h . 
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Table 1. Summary of the observing log and instrumental set up. 



Parameter 



13:00:11 13:00:02 12:59:54 12:59:45 12:59:36 12:59:28 
R.A. 



Fig. 1. Optical DSS image of the Coma cluster core. North is 
up, East is Left. The two red lines with the central gap represent 
the portion of the slit covered by the 2 CCDs in the FOCAS 
spectrograph. The scale in kpc given on the bottom left corner is 
computed assuming a distance of 97.7 Mpc, which is the average 
of the distances to NGC 4889 and NGC 4874 as given by NED. 



2. Spectroscopic observations and data reduction 

Long slit spectra were acquire d with the Faint Object Camera 
And Spectrograph (FOCAS, iKashikawa et all 120021) at the 
SUBARU telescope of the National Astronomical Observatory 
of Japan (NAOJ), on Mauna Kea (Hawaii, USA). Data were col- 
lected during two runs; in run #1 (April 2007) we obtained 8 
hours integration long slit data on the West side of the NGC 4889 
major axis, in the region in between NGC 4874 and NGC 4889, 
with a spectral resolution of 76 km s~' . In run #2 (May 2008) we 
obtained 5.5 hours integration of long slit data for the NGC 4889 
minor axis, with a spectral resolution of 96 km s _1 . In FigureQ] 
the observed field and slits locations are shown. In Table Q] we 
provide a summary of the observing log and instrumental set up. 
During these runs, we observed a set of kinematic template stars 
and Lick spectrophotometric standard stars for calibration to the 
Lick system. Long slit spectra on blank fields of the COMA 
cluster were also obtained for the sky background evaluation 
and to correct for large-scale illumination patterns due to slit vi- 
gnetting. In addition, we observed at least one flux standard star 
per night, to flux calibrate the spectra. 

Standard data reduction (bias subtraction, flat fielding to cor- 
rect for pixel to pixel chip sensitivity fluctuations) were per- 
formed using standard IRAF Q routines. Spectra distortion due 
to the optics were corrected using FOCASREE0 data reduction 
package within the IRAF environ ment. Cosmic rays were de- 
tected and removed using LACOS (van Dokkum 2001 ). Residual 
cosmic rays were removed by manually editing the spectrum. 
Wavelength calibration was performed using standard IRAF rou- 



2 IRAF is distributed by NOAO, which is operated by AURA Inc., 
under contract with the National Science Foundation. 

3 Information about the FOCASRED package can be obtained at the 
web page: http://www.naoj.org/Observing/Instruments/FOCAS/Detail/ 
UsersGuide/DataReduction/focasred.html 



Run 1 
13-14 Apr. 2007 



Run 2 
8-9 May. 2008 



Exp. time on galaxy 

Exp. time on sky 

Mean Seeing 

RA. slit center 

DEC. slit center 

Slit PA. 

Slit length 

Slit width 

CCDs 

Gain 

R.O.N. 

Pixel size 

Pixel binning 

Pixel scale 

Grism 

Filter 

Dispersion 

Observed range 

Calibrated range 

Calibration lamp 

Instrumental o" (,) 



16 x 30 min 
6 x 30 min 
~ 2" 
12:59:51 
+27:58:01 

81°7 
300 arcsec 
0.8 arcsec 
MIT 2.2K x 4K 
2.1 t I ADU 
4 e~ rms 
15 yum 
3 x 1 
ff.' 104 pixeF 1 (x3) 
VPH 450 
L600 
0.394 A pixeF 1 (xl) 
4200 - 5700 A 
4545 - 5665 A 
Th + Ar 
76 km s- 1 



1 1 x 30 min 
4 x 10 min 
~ 2" 
13:00:08 
+27:58:37 

-8°3 
300 arcsec 
1.0 arcsec 
MIT 2.2K x 4K 
2.1 e/ADU 
4 e~ rms 
15 /jm 
3 x 1 
tt' 104 pixeF 1 (x3) 
VPH 520 
L600 
0.400 A pixeF 1 (xl) 
4450 - 6050 A 
4700 - 5900 A 
Th + Ar 
96 km s-' 



Notes - w Computed at 5 100 A. 



tines. Comparison Th+Ar spectra were obtained during the night 
close in time to the scientific spectra. Because two different 
grisms were used for the two runs, we optimized the wavelength 
calibration in different spectral ranges: 4545A < A < 5665A, 
for run 1, and 4700A < A < 5900A, for run 2. With this 
choice, we were able to minimize the errors in the wavelength 
calibration, which are about +10 km s -1 as measured from the 
calibrated comparison spectra and from the sky line emission 
lines. Absorption lines as H/3, Mg I(AA5167, 5173, 5184) triplet 
and Fe 1{AA5270, 5328) fall inside the two selected wavelength 
ranges. Finally, the spectra were corrected for change in the slit 
illumination along the slit direction using the spectra obtained 
on sky blank fields. 



3. Sky subtraction 

Sky subtraction may introduce errors in the measurements from 
spectra of low surface brightness regions. To limit systematic 
effects as much as possible, we account for the sky subtraction 
in two ways. 

The first method is to obtain sky spectra from blank fields. 
The sky spectrum is then free of contamination either from the 
galaxy halos or intracluster light, but the spectra are not taken 
simultaneously with the scientific exposures. Therefore, if the 
relative intensities of the sky continuum, emission or absorption 
lines change with time, sky residuals are present in the galaxy 
spectra. 

The second method is to extract the sky spectrum from a 
region of the slit where the galaxy light is negligible (central part 
of the slit for run 1 , outer borders of the slit for run 2). The sky is 
now observed simultaneously to the galaxy observations but the 
disadvantage is that it might contain a residual contribution from 
the galaxy halo and/or intracluster light. We then may subtract 
spectral features, which belong to the galaxy and which we are, 
in fact, interested in measuring. 
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Fig. 2. Upper panel: difference between the mean sky spectrum 
measured at the offset position on the blank field and the sky 
spectrum measured at the long-slit center (run 1). The mean 
value of the difference is < 1 ADU. The positions where galaxy 
absorption line features should be visible (taking into account 
the galaxy radial velocity) are marked in green (see text for de- 
tails). The spectral ranges containing the most intense sky emis- 
sion lines are shaded in red. Lower panel: comparison of sky 
spectra obtained at different times during the night (black con- 
tinuous lines). The time variation of the relative intensity of sky 
emission and absorption lines is clearly visible. As in the up- 
per panel, the position of the galaxy Mg I absorption line is 
marked in green. For comparison, we show also the typical sky 
spectrum visible from the Kitt Peak National Observ atory {red 
dashe d line, rescaled to allow comparison, from lMassev & Folta 
2000). The same OH emission lines are visible. 



In the upper panel of Fig.[2]we compare the sky spectra ex- 
tracted using the two different methods. If the sky spectrum ex- 
tracted at the slit center (run 1) or at the slit edges (run 2) con- 
tains a small residual contribution from the galaxy (stellar con- 
tinuum and/or absorption line features) it should become visible 
when comparing with the sky spectrum extracted at the offset 
position. Residuals are approximately +0.5 counts (i.e. ~ 2% 
of the average sky level at 5100 A). This means that the level 
of the sky continuum computed with the two methods is the 
same, and differences between the two methods (if any) should 
be associated with spectral line features. However, we clearly 
see that larger residuals are present in the wavelength regions 
of the most intense sky lines. In particular, a residual spectral 
feature is seen around ~ 5280 A which is very close to the 
position of the Mg 1(^5167, 5173, 5184) triplet for the sys- 
temic velocity of NGC 4889. In the lower panel of Fig. [2] we 



compare the sky spectra extracted at offset position at different 
times during the night. We concentrate on the spectral region 
from 5150A to 5380A, which contains some OH emission lines 
(lOsterbrock et al.1 [T9 96). The variation of sky lines intensities 
over time is clearly visible and its amplitude (~ 1 — 2 AUDs) is 
consistent with the intensity of the spectral feature observed at 
~ 5280 A. These time variations indicate that the observed resid- 
ual at ~ 5280 A is related to the sky and can be best corrected 
for by using the sky observed simultaneously to the scientific 
observations. 



4. Extraction of spectra in low surface brightness 
regions 

Before measuring the stellar kinematics (Section|5]l and the line- 
strength indices (Section|6]l, we needed to bin spatially the spec- 
tra along the slit direction to increment the signal to noise ratio. 
Spectra from adjacent columns^] were added in order to reach a 
minimum signal to noise S /N ~ 10 for run 1, and S IN ~ 20 for 
run 2, which ensure reliable measurements of stellar kinematics 
and line strength indices. For run 1, the extracted spectra have an 
overall lower S /N than run 2 because the slit was sampling re- 
gions of lower surface brightness, and was illuminating regions 
of the CCD with bad cosmetics and hot columns which required 
masking. 

Summing adjacent spectra over part of the slit may affect the 
absorption line shapes by either smoothing the intrinsic velocity 
gradient or enlarging the absorption lines, if there is a substantial 
velocity gradient in the sampled region. 

Such an effect may be significant for run 1, where the slit 
sampled halo light of galaxies whose systemic velocities differ 
by ~ 700 km s -1 . We estimated the maximum number of adja- 
cent columns to be summed and extracted in a one-dimensional 
spectrum to be approximately ~ 30 - 40, depending also on the 
CCD cosmetics. Within this range, the maximum smoothing of 
the intrinsic velocity gradient is approximately ~ 20-30 km s -1 , 
which is the same order as the measurement errors. This value 
is obtained multiplying the number of columns by the maximum 
velocity gradient along the slit0, which we estimated to be ~ 0.8 
km s -1 pixel -1 . For run 2, we estimated the effect to be negligi- 
ble, because no velocity gradient was measured along the minor 
axis of NGC 4889 (Section©. 

In the case of NGC 4874, the ID spectra obtained from 40 
adjacent columns did not reach the required S IN. For those re- 
gions we used the following iterative procedure: 

1 . We selected N regions on the slit which are M < 30 pixels 
wide. Each region corresponds to a two-dimensional spec- 
trum, which we call "stripe" S,- (with i = \,...,N). Each 
of the N stripes is composed of M one-dimensional adja- 
cent spectra, with the wavelength direction oriented along 
columns. 

2. We extracted the one-dimensional spectrum S f, by summing 
the stripe along the slit direction, and assigned the luminos- 
ity weighted central position (in pixels) on the slit. We 
associated thus a spectrum S f and a central position Rf to 
each region of the slit. 



4 The spectra were oriented to have the dispersion direction along the 
vertical axis, and the slit aligned with the horizontal axis. 

5 This gradient was computed from the velocity difference measured 
at the slit edges (~ 800 km s -1 ) divided by the number of pixels in the 
slit (960). 
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Fig. 3. Examples of kinematic fits. Black: galaxy spectrum, 
green: portion of galaxy spectrum excluded in the fit, red: best fit 
template model. The S/N ratio per pixel of the spectra in panels 
from A to G are: 61, 15, 13, 10, 9, 8, 8, respectively. All spectra 
are normalized to the value of the best fit model at 5 100 A. 



3. We measured the radial velocity Vf of the spectrum S r. The 
radial velocity measurements are described in Section|5] 

4. We selected a reference spectrum S The other Af - 1 spec- 
tra were then shifted to the radial velocity VS of the refer- 
ence spectrum. The AV applied to the S. ,, spectrum was 
AV = 6V ■ (Rf - R%, f ), where 6V (in km s -1 pixel -1 ) is the 



NGC 4889 major axis R = 65 kpc 
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Fig. 3. - Continued. 

(unknown) velocity gradient along the slit. The shifted spec- 
tra were then combined together to form the total spectrum 
gTOT j n tne ra( jj a j b m centered at R c ... 6V was determined 

ref 

via an iterative procedure that minimized the line widths in 
the combined spectrum S T0T . 

The above scheme and the iterative procedure for the deter- 
mination of 6V assumed a constant velocity gradient along the 
slit. In order to minimize systematics that may arise from this 
assumption, we limited the radial binning to a maximum of 3 
adjacent stripes 5Vi,;,i+i, each of them 30 pixels wide at most. 

5. Stellar kinematics 

Stellar kinematics from the extracted spectra in the radial 
bins w ere measured using the "Penali zed Pixel-Fitting method" 
(PPXF, ICappellari & Emsellernl l2004h. A library of kinematic 
template stars (MILES library. ISanchez-Blaza uezet al.l [2006) 
was used together with stellar templates observed during the ob- 
servations. The template spectra of the MILES library were con- 
volved with a Gaussian function to match the FOCAS spectral 
resolution. We show some of the fit results in Figure [3j the in- 
nermost and outermost spectra extracted along the minor axis of 
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Fig. 4. Radial velocity {upper panels) and velocity dispersion 
profiles (lower panels) along the slit in run 1 (velocities are 
expressed in the heliocentric reference system). Measurements 
along the slit are divided in two panels, with left panels refer- 
ring to the part of the slit close to NGC 4889 and right pan- 
els referri ng to the part of the slit close to NGC 4874. Crosses: 
data from Me hlert et al.l (2000) shi fted to the sy stemic velocities, 
6495 km s 1 for NGC 4889 dMoore et al.l2002l) and 7205 km s 1 
for NGC 4874 (ISmith et al.ll2004l) . Filled circles: measurements 
from this paper for run 1 . 
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Fig. 5. Stellar kine matics along the mi nor axis of NGC 4889. 
Crosses: data from lCorsini et al.l (12008). filled circles: measure- 
ments from this paper for run 2. 



NGC 4889, the three outermost spectra along the major axis of 
NGC 4889 and the two spectra of NGC 4874 are also shown. 

The wavelength ranges used for the fitting are 4600A < A < 
5550A for run 1, and 4800A < A < 5850A for run 2. Bad pixels 
coming from residual cosmic rays or emission lines were prop- 
erly masked and not included in the fitting procedure. 

For each radial bin, the PPXF code built the optimal template 
which represents the observed spectrum by a linear combination 
of different stellar templates. The use of a large spectral library 
and the creation of one optimal template for each radial bin min- 
imized errors in the kinematics due to template mismatch. In the 
outermost bins (R > 60 kpc in run 1 and \R\ > 20 kpc in run 2) 
the construction of the optimal stellar template in PPXF failed, 
because the absorption line features are less pronounced. The 
kinematics in these regions were measured using the best tem- 
plate obtained from the nearest bin in which it could be otained. 
Referring to Fig. [3] spectra in panels D, E, F and G were fitted 
using the optimal template determined when fitting the spectrum 
in panel C. The spectrum in panel B was fitted with the optimal 
template obtained fitting the spectrum at 16 kpc (not shown in 
the figure). This approximation might introduce some system- 
atic errors in the kinematics of the outer bins due to template 
mismatch. Nevertheless, in the radial bins (\R\ < 20 kpc in run 
2), we did not observe significant differences in the kinemat- 
ics when either the "central" optimal template at R = -4.9 kpc 
or the "local" optimal templates, found by PPXF for those bins, 
were used. This confirmed that errors on the kinematics caused 
by template mismatch were negligible for the inner bins, and 
gave us confidence that the template mismatch is not important 
in the outer bins. Errors on the measured kinematics were deter- 
mined by means of Monte Carlo simulations, analyzing spectra 
generated from the best fit template model and adding the appro- 
priate noise. 

The measured radial velocities were then shifted to the helio- 
centric reference system using the MIDAS^task COMPUT/BARY. 

The fitted me an velocity, velocity dispersion and Gauss- 
Hermite moments (Gerhard 1993; van der Marel & Franx 119931) 
and their errors are shown in Figures |4] and [5] for runs 1 and 2, 
respectively, and are listed in Table [2] Gauss-Hermite moments 
were fitted only for spectra within \R\ < 20 kpc along the minor 
axis. 



6. Line strength indices 

Line strength indices were measured for NGC 4889 using 
observations from both runs. The spectra of NGC 4874 did 
not have enough S /N for this. We measured the H/3, Mgi, 
M g2, Mg b, Fe 5 270 an d Fe 5335 line-strength ind ices as defined 
by iFaber et al.l (1 19851) and IWorthev et alt (11994 . The spectra 
were convolved with a Gaussian function to match the spec- 
tral _£e^o2ution_£f_theLic^_sy stem (FWHM = 8.4 A at 5100 
A. IWorthev & Ottavianil Il997h . Three Lick standard stars ob- 
served during the runs were used to correct for offsets to the 
Lick system. The measured offsets (FOCAS - Lick) are: 0.03 A 
(H/3), -0.30 A (Fe5105), -0.028 mag (M gl ), -0.016 mag (Mg 2 ), 
-0.037 A (Mg b), -0.36 A (Fe5270), 0.024 A (Fe5335) and 0.024 
A (Fe5406). 

No focus correction was applied because atmospheric seeing 
was the dominant effect during observations. 



6 MIDAS is developed and maintained by the European Southern 
Observatory. 
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Fig. 6. Li ne-strength Lick ind ices for NGC 4889 me asured along the major (left panel) and the minor axis (right panel). Crosses: 
data from Meh lert et alj |2000 (major axis) and lCorsini et alj|2008l (minor axis). Filled circles: measurements from this paper. Solid 
lines indicate the systematic deviations of the line-strength measurements caused by adding different residual sky levels to the 
extracted galaxy spectrum. From top to bottom, the residual sky levels we added are +2% (i.e. sky subtraction is underestimated by 
2%), +1%, +0.5%, -0.5%, -1%, -2% (i.e. sky subtraction overestimated by 2%). See Section[3]for details. 



Errors on line-strength indices were determined via Monte 
Carlo simulations, which accounted for the errors on radial ve- 
locity measurement also. Errors from Monte Carlo simulations 
are consisten t with those predicted by the analytic formulae of 
ICardiel et all (fl998r) . 

In the following we describe the details of additional cor- 
rections we applied to the measurements in order to take into 
account instrumental effects. The final measured values are 
liste d in Table [2] and a re shown in Figure [6] together with 



the Mehlert et al. (2000) data along the major axis, and the 



ICorsini et al. (120081) data along the minor axifl 



6. 1 . Correction for intrinsic broadening 

The line-strength measurements were corrected for the galaxy 
intrinsic broadening fo llowing the procedure described in 
iKuntschner et alJ d2006l) . Coefficients C; for each index / were 
determined as the ratio of i) the "intrinsic" values (Iq, in 
Angstrom) measured on the optimal stellar template; and ii) the 
line-of-sight-velocity-distribution "modified" values (/losvd) 
measured on the best fit model (i.e. the optimal template con- 
volved with the galaxy LOSVD). Corrected galaxy line strength 
indices (in Angstrom) were then given by: 



= C, ■ I n 



^LOSVD 



(1) 



7 We show the mean iron < Fe >= | (Fe5270 + Fe5335) radial pro- 
files and not the profiles of both iron indices to facilitate the comparison 
with the literature, where only < Fe > is provided. 



6.2. Correction for residual emission lines 

The long-slit spectra for the NGC 4889 minor axis showed 
strong unresolved emission lines in the wavelength regions ad- 
jacent to the Mgl bandpass and the H/3 blue pseudo-continuum. 
These strong lines were at wavelengths that are not associated 
with known galaxy emission lines or sky/auroral lines. Our con- 
clusion was that they may be associated with internal reflec- 
tions in the spectrograph. To measure the spectral indices in 
the regions affected by the unidentified emission lines, we si- 
multaneously fit the galaxy spectrum with a stellar template 
plus Gaussian functions at the wavelength of the detected emis- 
sions. The best-fit Gaussians were then subtracted from the spec- 
trum, to remove the observed emission lines. As a further test, 
we checked the measured values of the line-strength indices 
against the measurements on the best-fit stellar template, con- 
volved with the galaxy line of sight velocity distribution deter- 
mined in Secfion[5] The two sets of values were consistent. 

6.3. Systematic effects caused by sky subtraction 

We then quantified the systematic errors in the line strength in- 
dices measurements caused by sky subtraction . A residual sky 
contribution of +0.5%, +1% and +2% of the subtracted sky 
spectrum^] was added to the extracted galaxy spectra, and then 
we performed measurements of the line-strength indices on these 
modified spectra. Within the errorbars, these new values were 
consistent with the measurements performed on the best back- 



8 Residual sky contribution larger than 2% would produce detectable 
sky features either in emission or absorption in the galaxy spectra. 2% 
is also the maximum difference between the continuum of the sky mea- 
sured at blank fields and the one measured simultaneously to galaxy 
observations (see Section[3]and Fig.f2]for discussion). 
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ground subtracted galaxy spectra, but for the Mg2 index. This 
was probably due to the fact that the Mg2 pseudo-continua were 
defined on spectral regions not adjacent to the index passband 
(unlike for the other spectral indices), therefore variations of the 
sky level have a larger effect on this index. 

7. Results and discussion 

In NGC 4874 the measured stellar kinematics extend out to ~ 50 
kpc (~ 1.4/? e ) along its East side. Measurements of line strength 
indices were not possible. 

In NGC 4889 our measurements extend out to 65 kpc (~ 
43R e ) along its West side (close to its major axis) and ~ 35 kpc 
(~ 2.3R e ) along its photometric minor axis. This represents the 
most spatially extended datasets with both stellar kinematics and 
line strength indices for a brightest cluster galaxy. 

7.1. Kinematics 



2009; 

iProctor et alJ l2009t ISchuberth et alJ I20091 McNeil et al. 2010 
A&A submitted). These galaxies do not show an increase in the 
velocit y dispersion profile as observed in other BCGs (e.g NGC 
6166 in lKelson etatll2002h . which was interpreted as generated 
by stars not bound to the galaxy but free floating in the cluster 
potential. 

We can thus compare the kinematic properties of the outer 
regions of NGC 4874 and NGC 4889 with the kinematic proper- 
ties of the outer halos of other early-type galaxies. To do this, we 
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1.2 7 

< 0.8 r 
V 0.6 7 
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0.2 7 



42 7 
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NGC 4874. In NGC 4874 the outer mean velocity along the line 
of sight is < V >= 7150 ± 37 km s _1 , about ~ 55 km s _1 lower 
than the gal axy systemic velo city (V syst = 7205 + 19 km s _1 , as 
reported by Smith et al. 2004). It is not straightforward to asso- 
ciate it with rotation or tidal effects, given the errors on the cur- 
rent measurements and the lar ge uncertainties on the galax y sys- m 
temic velocity. As examples, Ide Va ucouleur setall (1 1 99 it RC3 
hereafter) report a systemic velocity of 7152 + 16 km s _1 , while 
NED reports a systemic velocity of 7224 +11 km s _1 . The mean 
outer velocity dispersion in NGC 4874 is cr = 283 + 33 km s' 1 , 
which is consistent with the central value found by Smit h et al.l 
(2004), suggesting a flat velocity dispersion profile with radius. 
NGC 4889. In NGC 4889 the mean outer velocity is < V >= 
6326 + 44 km s _1 , which is ~ 169 km s _1 lower than the galaxy 
systemic velocity (V CTJ> = 6495 + 13 km s _1 , as reported by 
Moor e et"a"D 12002). The measured velocity offset is significant 
when compared with measurements errors and uncertainties on 
the systemic velocit y. Although rota t ion in the halos of e llipti- 
cals is not rare (e.g. lRix etal1ll999t ICoccato et"al1 12009). it is 
not straightforward to consider the velocity offset as evidence 
of rotation, since we miss measurements on the East side of the 
galaxy. The outer mean velocity dispersion we measured in NGC 
4889 is cr = 266 + 17 k m s' 1 , consis t ent w ith the outermost 
measurements (~ 10") by iMoore et alJ (120021) . revealing a flat 
velo city dispersion profi le outside the central drop. 

iGerhard et al.l d2007l) developed a model for the core of the 
Coma cluster whereby the two BCGs, NGC 4874 and NGC 
4889, are in the process of merging, tidally stripping each oth- 
ers' halos to create an elongated distribution of intracluster light 
in which they are embedded. It is worth noting that, out to the 
radii probed by the present observations, the velocity dispersion 
profiles of NGC 4874 and NGC 4889 are nearly constant and in- 
dicative of stars still bound to the central galaxy. This is consis- 
tent with recent velocity dispersion profiles measurements based 
on planetary nebulae and globular clusters for local brightest 
cluster galaxies (e. g. Coccato et al. 120091: iDohertv et alJ l2( 
r et alJ 120091: ISchuberth et al.l 120091: McNeil i 
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Fig. 7. Correlations between the outermost measured value of 
the halo velocity dispersion and mean < V/cr > of the halo 
(upper panel), total X-ray luminosity (central panel) and total 
luminosity in the B band (lower panel). Crosses are data from 
ICoccato et all |2009), filled circles are data from this paper. 



refer to the results of lCoccato et al.l (120091) based on radial velo- 
cities of Planetary Nebulae in the halos of early-type galaxies. 
Their sample covers the range 50 km s _1 < cr < 220 km s _1 of 
velocity dispersion in the outer halo. The two galaxies studied 
in this pa per fall in the upper cr range of the galaxy sample ex- 
plored bv ICoccato et al.ll2009l and t his e nables us to extend their 
cr range. Following ICoccato et al. (2009), in Figure Q we com- 
pare the halo velocity dispersion with other galaxy properties, 
< V/cr > (we assume that the observed velocity offset is due to 
rotation), total X-ray luminosity (Lx = 6.8- 10 41 erg s~ | for NGC 
4874 and L x = 5.9 ■ 10 42 erg s 1 for NGC 4889, from lPeUegrinil 
120051) and total fi-band magnitude (B T = -22.45 for NGC 4874 
and B T = -22.50 for NGC 4889, from RC3). Both NGC 4874 
and NGC 4889 follow the same trend found in other early-type 
galaxies, therefore supporting the earlier statement that the stars 
in their halos are still bound to the inner galaxies in the radial 
range explored in this work. 

7.2. Line strength indices 

In NGC 4889, the measured set of Lick indices provides infor- 
mation in the radial interval from ~ 7 to ~ 65 kpc, and is comple- 
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mentary to data available for the central regions dMehlert et al.1 
l2000HCorsini et al.ll2008h . 

In the region of overlap the current data set and the published 
values agree along the minor axis (\R\ ~ 5 kpc), but for H/3. 

In the central p arts of the galaxy, the mean velue from 
ICorsini et al.l (|2008) data is higher than the inward extrapolation 
of our data by 0.3 A. Their outermost central value differs from 
ours of 1.1 A at the same radius (R ~ -5 kpc, see Pig. [6j>. 

The difference could be related either to the limited number 
of stars we used to correct the FOCAS observations to the Lick 
system, or to sky subtraction, affecting Corsini's measurements 
(see Section |6}. But in these cases we would expect a discrep- 
ancy also in other indices, and not only in H/3. Another reason 
could be the presence of emission lines, but neither in our spec- 
tra nor in Corsini's spectra H/3 emission lines were detected at 5 
kpc. Therefore, the reason for the discrepancy remains unclear. 

Values of H/3 range from 0.5 to 2.3 A. Line strength indices 
related to Mg I and < Fe > show declining profiles, from the 
central values of Mg 2 ~ 0.37 A, Mg b ~ 5.5 A and < Fe > ~ 3.4 
A, to outer values of Mg 2 ~ 0.15 A, Mg b ~ 3 A and < Fe > 
~ 2.5 A. 

In lCoccato et al.l J2010) the combined set of line strength in- 
dices for NGC 4889 is used to derive age, metallicity and alpha- 
enhancement for single stellar population models and their radial 
gradients from the center out to 60 kpc. There we find evidence 
for a change in the stellar population gradient and content, sug- 
gesting a different evolutionary history between the galaxy and 
its halo. We refer to this forthcoming paper for a detailed discus- 
sion of the implications of these measurements for the properties 
of the stellar populations in NGC 4889, and for the constraints 
they put on the galaxy formation mechanism. 
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Table 2. Stellar kinematics and line strenght indices for NGC 4874 (run 1) and NGC 4889 (runs 1 and 2). 



R V e a h 3 h 4 Hj8 Fe5105 M gl Mg 2 Mg b Fe5270 Fe5335 Fe5406 

(kpc) (kms- 1 ) (kms- 1 ) (A) (A) (mag) (mag) (A) (A) (A) (A) 

0) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
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6315 ± 18 


244 
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2.91 ± 0.39 


2.43 
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0.38 


2.71 
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0.46 






- 




60.31 


6348 ± 24 
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± 24 
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1.96 
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0.34 
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0.02 ± 0.09 


0.144 ±0.010 


2.65 ± 0.40 


2.12 
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0.46 






- 
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6323 ± 27 
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± 28 
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2.08 w = 


.0.33 
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0.03 ± 0.09 


0.186 ±0.011 


3.93 ± 0.40 
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6447 ±8 


310 ±9 


0.02 ± 0.28 


0.00 ± 0.26 


1.39 


± 


0.16 


4.50 ± 0.46 


0.18 ±0.04 


0.285 ± 0.005 


4.43 ±0.19 


2.34 


± 
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2.34 


± 


0.18 


1 


39 


± 


0.21 


-7.57 


6449 ± 13 
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± 14 


0.04 ±0.12 


-0.03 ±0.11 


1.36 


± 


0.31 


4.53 ± 0.33 


0.18 ±0.08 


0.287 ± 0.010 


4.43 ± 0.36 
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± 
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± 


0.35 


1 


44 


± 


0.15 


-4.92 


6447 ± 9 
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Notes - Measurements marked by w were obtained by removing an emission line component from the spectrum. The symbol - means that the measurement has not been performed. 



